Real-time renormalization group and charge fluctuations in quantum dots 
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We develop a perturbative renormalization-group method in real time to describe nonequilibrium 
properties of discrete quantum systems coupled linearly to an environment. We include energy 
broadening and dissipation and develop a cutoff-independent formalism. We present quantitatively 
reliable results for the linear and nonlinear conductance in the mixed-valence and empty-orbital 
regime of the nonequilibrium Anderson impurity model with finite on-site Coulomb repulsion. 

72.10.Bg, 73.23.Hk, 05.10.Cc, 05.60.Gg 



CD 



I 

O 

o 



> 

O 

oo 
O 



I 

o 
o 



Renormalization group (RG) methods are standard 
tools to describe various aspects of condensed-matter 
problems beyond perturbation theory . Many impurity 
problems have been treated by numerical RG with excel- 
lent results both for thermodynamic quantities and spec- 
tral densities [^J^. These RG techniques, however, can- 
not describe nonequilibrium properties like the nonlinear 
conductance, the nonequilibrium stationary state, or the 
full time development of an initially out-of-equilibrium 
system. To address these aspects we present here a per- 
turbative RG method, formulated for strongly-correlated 
quantum systems with a finite number of states coupled 
linearly to external heat or particle reservoirs. Funda- 
mentally new, we work on a Keldysh contour and gen- 
erate non-Hamiltonian dynamics during RG, which cap- 
tures the physics of finite life times and dissipation. Fur- 
thermore, no initial or final cutoff in energy or time space 
is needed, i.e., large and small energy scales are accounted 
for correctly like in flow-equation methods |4| . Although 
correlation functions can also be studied, physical quan- 
tities like spin and charge susceptibilities or the current 
can be calculated directly without the need of nonequi- 
librium Green's functions. 

The purpose of our RG technique is to describe quan- 
tum fluctuations which are induced by strong coupling 
between a small quantum system and an environment. 
There are several recent experiments which show the 
importance of quantum fluctuations in metallic single- 
electron transistors |^ and semiconductor quantum dots 
|^-|| (see ||] for an overview over theoretical papers). 
Due to the renormalization of resistance and local en- 
ergy excitations, anomalous line shapes of the conduc- 
tance have been observed, which can not be explained 
by golden-rule theories. Therefore, the RG approach pre- 
sented here can treat simultaneously strong coupling to 
the reservoirs, strong Coulomb interaction on the island, 
and finite bias voltage. We apply the technique to a 
quantum dot with one spin-degenerate state and present 
for the first time quantitatively reliable results for the 
nonlinear conductance in the mixed-valence and empty- 
orbital regime. For applications to the single-electron box 



and the single-electron transistor we refer to Ref. (lOj . 

We consider the nonequilibrium Anderson impurity 
model with finite on-site Coulomb repulsion U. It con- 
sists of a single, spin-degenerate state in a quantum dot 
which is coupled via tunneling to a left and right reser- 
voir. The Hamiltonian H = Hj.es + Hd -\- Hx consists 
of three parts. Hj-es = Y.kar <^kara\„r'^kar dcscribcs two 
reservoirs, r = L, R, where k labels the states and a de- 
notes the spin. The isolated quantum dot is given by 
Hn = J2a^<^^<^ + Un^ni with 11^ = c\ca- The coupling 
to the reservoirs is described by the standard tunneling 
Hamiltonian Ht = Sfeo-rC^fc'^Lr'^CT + h.c). The goal is 
to evaluate the reduced density matrix of the dot sys- 
tem, poit) = TrresPit), and the expectation value of the 



current operator I = II = iej^kai'^t 
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h.c). 



We start with the time evolution of the dot system, 
pjy{t) = Tr^es {e~'-^V(0)e*^*}- Initially, the system 
is decoupled, and each reservoir is in thermal equilib- 
rium, p(0) = pd{0)p'i^ p'r ■ Nonequilibrium situations 
arise when the two reservoirs have different electrochem- 
ical potentials /i^. As in Refs. |^Jll[l we expand the for- 
ward/backward propagators exp(=pii7i) in Ht, and per- 
form the trace Tr^es by applying Wick's theorem with 
respect to the reservoir field operators. All terms can be 
represented diagrammatically as shown in Fig. |l|. Tun- 
neling vertices are ordered along a closed Keldysh con- 
tour. They are connected in pairs by the contractions 
(dashed lines in Fig. 0) j+it) = |T,r|2(aL.(i)«fe..> 



or 7,, (t) = J2k l^fe I (afe<Tr(<)aLr), depending on the rel- 
ative time-ordering of the vertices on the contour. We 
obtain for r = L, R, 
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2/3sinh[7r(t-iO+)//3] 



(1) 



with 13 = llikeT), and F, = 2^Y.k \TL?KE - ^kar)- 
The solid line in Fig. ^ represents free time evolution of 
the dot. As a result, we have obtained an effective theory 
of the dot degrees of freedom while the reservoirs have 
been integrated out. 

To derive a kinetic equation we will call diagrams ir- 
reducible if any vertical cut crosses at least one dashed 
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line. Fig. ^shows four such irreducible blocks (a)-(d). We 
denote the sum over all irreducible diagrams between t' 
and t by the kernel T,{t — t'). The full time evolution, 
i.e., a sequences of E blocks, can be described by a self- 
consistent equation. Differentiating with respect to time 
t leads to the standard kinetic equation [^11 



PD{t)+lLDPD(t)= f dt'^{t-t')pD{t') 

Jo 
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and E are superoperators acting on 

are 



Here, Ld = [Hd, 

PD{t), i.e., the 4** = 256 matrix elements E^^^'^ ^^s: 
labeled by four dot-states with si/2 {s'if2) referring to 
the forward (backward) propagator. The l.h.s of Eq. (|^) 
describes the time evolution of the dot in the absence of 
tunneling whereas the r.h.s contains the dissipative part 
which drives the dot distribution into a stationary state. 
In Laplace space Eq. (^) reads [z — Ld — iE(z)]/5£)(z) = 
«P£)(0). Thus, the knowledge of E(z) provides the full 
time evolution of the dot distribution. The stationary 
solution follows from [Ld + iE(0)]p|5 = 0- The differ- 
ent electrochemical potentials of the reservoirs (leading 
to a nonequilibrium stationary state) enter via the pair 
contractions (|l|) in the kernel E(z). 

Sequential tunneling or the golden rule results fl^ j can 
be recovered by calculating the kernel E(z) perturba- 
tively in first order in the tunneling coupling . The aim 
of the present paper, however, is to go beyond and cal- 
culate the kernel nonperturbatively by a systematic RG 
procedure. The idea is to integrate out all contraction 
lines one after another, beginning with the shortest one. 
Fig. shows four examples (a)-(d) in second order in F. 
In these examples, the lines with index R are integrated 
out first. This leads to a renormalization of the dot prop- 
agator in (a) and the tunneling vertex in (c). In (b) and 
(d) the contraction line R connects the forward with the 
backward propagator. This element does not occur in any 
equilibrium theory but arises for a Keldysh contour in a 
natural way. Such contractions lead to non-Hamiltonian 
dynamics for the dot distribution and, therefore, account 
for dissipation. To treat these cases in the same way as 
(a) and (c), it is convenient to view the forward and back- 
ward propagator formally as one double line, see Fig. 
As a consequence, (b) and (d) fall into the same topolog- 
ical classes as (a) and (c), respectively. The price is that 
a "state" on the double line has to be specified by two 
dot states, one for the upper and one for the lower prop- 
agator. This leads to the superoperator matrix notation 
used in Eq. ^ . The tunneling vertices on the double 
line act on these double states. There are 16 different 
vertices, specified by four indices, p and p = rjar, where 
p = ± indicates whether the unrenormalized vertex acts 
on the forward/backward propagator, 77 = ± represents 
tunneling in/out, a i is the spin, and r — L, R refers 
to the reservoir index. As shown in example (d) of Fig. 
after being renormalized a vertex will in general act on 



both the forward and backward propagator simultane- 
ously. Double lines without vertices corresponds to a free 
time evolution exp(— iL/jt) of the dot (with renormalized 
Ld)- Finally, we denote the rightmost (leftmost) vertex 
of the kernel E by {Rp). They are renormalized in a 
different way than C^. 

The flow parameter tc in our RG procedure is the 
largest length of those contractions that have already 
been integrated out. An infinitesimal RG step is estab- 
lished by integrating out all contractions with a length t 
between tc and tc + dtc- For tc 00 all contractions are 
integrated out. We find the RG equations (summation 
over repeated indices implicitly assumed) 
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-lZ{tc)I,,:Al{tc)Bl, , 

-^7r.'(^c)vc^(ic)c^; , 
Cdt [pp'j,,,^ci\it)ci 
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r ' {tc) 



dt 



Pp[Ip,p'CP\{t)CPB^^it-tc) , (7) 



We introduced 7^^,(i) 



i^^.y-j Ir^ '^ip't)Sp,ti' for the contrac- 
tion between two tunneling vertices on the double line 
(/i = rjar, p = —rjar, p' ^ 
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= rj'a'r') to get a com- 
The interaction picture is defined by 
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B^e - The diagonal superoperator /j 
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pact notation. 
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together with the factors pp'i account for possible mi- 
nus signs due to Fermi statistics. We get {Ipp')ss' .ss' — 
{pp')^"^^"' , where Ng is the particle number for state 
s- All terms arise naturally from our intuitive picture 
set up in Fig. ^ except for the second term on the r.h.s 
of Eqs. (||) and (||). The latter are correction terms to 
account for correct time ordering fl^ . 

To evaluate the average current we repeat the above 
derivation of the RG equations but replace the un- 
renormalized boundary vertex operator A^^ by the cur- 
rent vertex /. The analog of Eq. (||) gives the cur- 
rent kernel E/(z) and the current follows from {I){z) — 
Ttd'^i{z)p{z)- 

Taking matrix elements, all integrals in Eqs. (^)-(0) 
can be calculated analytically, and pure differential equa- 
tions are left which we solve numerically. Most of the ma- 
trix elements which are initially zero remain unchanged 
under the RG flow. We end up with 26 equations for E 
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and Lo each, and 192 equations for C, A, and B each to 
solve. 

The RG flow of conventional poor man scaling and op- 
erator product expansion methods stops at some " charac- 
teristic" time (or inverse energy) scale. All results, there- 
fore, depend on a low-energy cutoff. This is not the case 
in our formulation since we have not expanded the propa- 
gation exp(±jL£)tc) in tc (it occurs in the interaction pic- 
ture of the vertex operators) and can, therefore, integrate 
out all time scales. Furthermore, we work on a Keldysh 
contour and find that the renormalized Lu-operator of 
the dot system can no longer be represented as the com- 
mutator with a renormalized Hamiltonian Ljj ^ [Hd, ■]. 
This corresponds to non-Hamiltonian dynamics and is 
needed to describe the physics of dissipation and finite 
life times. 

An important question concerns the validity range of 
the present approach. Although no general statement is 
possible, the RG flow itself will indicate the breakdown of 
the procedure by instabilities or increasing coupling con- 
stants. Due to the neglect of higher-order vertex correc- 
tions, one might not trust the procedure for too large cou- 
pling constants. However, for the specific example of the 
Anderson model in the mixed-valence regime (see below) 
we reproduce exact Bethe ansatz solutions within high 
accuracy even in the strong-coupling regime. The same 
was achieved for the ground-state energy of the single- 
electron box |10|. A possible explanation is the gener- 



ation of rates and energy broadening during RG which 
improves the accuracy and sets another energy scale for 
the cutoff of the RG flow (for the present problem given 
by r) . This is consistent with our numerical solution and 
with a recent analysis of poor man scaling on a Keldysh 
contour in Ref. [Q. 

In Fig. ^ we show the linear conductance G for the 
Anderson model as function of e — e^, which is var- 
ied experimentally by the gate voltage. For T = and 
e > — 0.4r, the deviation from the exact Friedel sum 
rule G\t=o = {2e^ / h) sm^ {-k {n) / 2) is less than 2%. Here, 
r/2 = Tl — Ffl and (n) = ^^^("'cr)- ^^'-^ same regime, 
the average occupation (ri) agrees with the exact Bethe 
ansatz solution [|l5| within 3% (inset of Fig. ^). Since 
the RG flow is strongest if both temperature and bias 
voltage vanish, this comparism gives good support that 
our results are reliable in the whole mixed-valence and 
empty-orbital regime for all temperatures and bias volt- 
ages. Furthermore, we find excellent agreement with the 
exact Bethe ansatz solution for the magnetic and charge 
susceptibility at e = as function of T and J7, respec- 
tively . Thus, the neglect of higher-order vertex terms 
seems to be justified. 

In the Kondo regime, i.e., for e < — 0.4r, spin fluc- 
tuations arc dominant and double-vertex terms become 
important. In this regime, deviations from the Friedel 
sum rule can be seen in Fig. ^ Experimentally, such de- 
viations occur as well [l6|-EI, since it is difficult to reach 



the T = limit. Our theory accounts for charge fluctua- 
tions very well which are responsible for the temperature- 
dependent renormalization of the dot level as shown by 
the shift of the maximum peak in Fig. ^. The fact that 
the conductance increases along the zero temperature re- 
sult on the r.h.s of the conductance peak is in good agree- 
ment with experiments (6|Q|. A fit of the temperature 
dependence of the linear conductance to the experiment 
in Ref. |^ is shown in Fig. ^ for various values of e. The 
agreement is quite well except for T ^ O.SF where the 
experimental conductance is higher. A reason for this 
discrepancy might be the presence of other levels with a 
level spacing of order F. For a low- lying level the conduc- 
tance increases monotonically with decreasing tempera- 
ture and finally saturates. The saturated value at T = 
is not identical to the unitary limit {e^ / h)&^ [X n/T^ since 
we are still in the mixed-valence regime. However, the 
value agrees perfectly with the Friedel sum rule as al- 
ready stated above. In the empty-orbital regime, e ^ 0, 
the conductance shows a local maximum by varying T. 
Here, the renormalized e is above the Fermi level and will 
be occupied at high enough temperature. This increases 
the conductance with increasing T. 

Fig. |showsG = d//(iX^ as a function of the bias volt- 
age for e = — 0.4r. We obtain a zero-bias maximum with 
an amplitude increasing monotonically with decreasing T 
(see Fig. ^ . This maximum is not due to the presence of 
a Kondo resonance but stems from the renormalized dot 
level being near the Fermi levels of the reservoirs. The 
width of the peak in Fig. |5| scales with F which is the 
relevant energy scale in the mixed-valence regime. Since 
the zero-bias maximum has a symmetric line shape, our 
results can be used for a fit of the important energy scale 
r in experiments. 

In summary we have presented a new RG technique to 
calculate nonequilibrium properties of dissipative quan- 
tum systems beyond perturbation theory. The method is 
very flexible and was applied to the Anderson-impurity 
model with finite U . In the mixed-valence and empty- 
orbital regime the Friedel sum rule and the exact Bethe 
ansatz solution were reproduced at zero temperature. 
The results are consistent with experiments. 
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FIG. 3. Linear conductance as a function of the dot level. 
T/r = 0,0.05,0.1,0.25,0.5 (solid lines from top to bottom) 
and U = 6r. The inset shows the average occupation for 
T = 0. 
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FIG. 4. Linear conductance as a function of temperature 
in comparism to experiment i7 = 6.2r, SFlFr/F^ = 0.6, 
and r = 3423mK. 
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FIG. 1. Example of a diagram on the Keldysh contour for 
the matrix element pD(t)oo- The contractions are labeled by 
the reservoir indices L and R. 
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FIG. 2. The same figure as Fig. hi but the two lines taken 
together. 
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FIG. 5. Nonlinear conductance as a function of the bias 
voltage. 
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